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ABSTRACT 


We  are  interested  in  classifying  the  solutions  of  Riemann  problems  fbr 
the  2*2  conservation  laws  which  have  homogeneous  quadratic  flux 
functions.  Such  flux  functions  approximate  an  arbitrary  2*2  system  in  a 
neighborhood  of  an  isolated  point  where  strict  hyperbolicity  fails.  This 
problem  was  motivated  by  Marchesin  and  Paes-Leme  who  discovered  such  a 
singularity  in  a  system  of  equations  arising  in  oil  reservoir  simulation  , 
[15,17].  In  [18],  Schaeffer,  Shearer,  Marchesin  and  Paes-Leme  solved  the 
Riemann  problem  for  this  system  in  a  neighborhood  of  the  singular  point.  In 
[8] ,  Isaacson  and  Temple  outlined  a  program  for  classifying  such  singularities 
by  means  of  locating  normal  forms  for  the  equivalence  classes  of  equations 
generated  by  linear  changes  in  dependent  variables.  A  2-parameter  family  of 
such  normal  forms  was  found  by  Plohr.  In  the  important  work  of  Schaeffer  and 
Shearer  [17],  a  new  norr.^l  form  was  found  which  reduced  the  classification  of 
integral  curves  to  a  theorem  of  Darboux  on  the  classification  of  umbilic 
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points  for  honogeneous  cubic  equations^  %a  Integral  curves  fall  Into 

four  Isooorphlsm  classes,  called  Regions  I-IV.  In  this  paper  we  give  the 
solution  of  the  Rlemann  problem  for  the  systems  In  Region  XV  which  exhibit  up- 
down  syimnetry.  A  presentation  of  the  solutions  of  the  corresponding  systems  In 
Regions  II  and  III  will  follow. 

Our  analysis  uses  a  numerical  determination  of  the  Hugonlot  loci. 
Moreover,  a  new  phenomenon  occurs  In  Regions  11-IV  that  does  not  occur  In 
Region  I:  lines  exist  on  %«hlch  one  of  the  eigenvalues  has  the  value  of  the 
eigenvalue  at  the  singular  point.  Bifurcations  occur  at  these  lines,  and  the 
dyn^unlcs  of  the  lines  as  the  parameters  In  the  normal  form  are  varied  gives  a 
geometric  Interpretation  of  the  boundaries  between  Regions  I-III.  The  program 
for  classifying  hyperbolic  singular  points  In  2x2  systeais  of  conservation 
laws  Is  being  carried  out  jointly  by  the  authors  named  above. 

AMS  (MOS)  Subject  Classifications:  65M10,  76N99,  35L65,  351.67 

Key  Words;  Rlemann  problem,  non-strlctly  hyperbolic  conservation  laws, 
umblllc  points 

Work  Unit  Numbers  1  (Applied  Analysis)  and  3  (Numerical  Analysis  and 

Scientific  Computing) 


SIGHIFICAMCE  AMD  EXPIAMATIOM 


A  2x2  system  of  conservation  laws  is  a  system  of  partial  differential 
equations  of  the  form 

(1)  +  f(u)  *  0 

where  u  ~  (U/V),  f  >  {f,q)  are  in  R  x  k  and  where  x  £  R,  t  >  0.  Such 

equations  arise  in  gas  dynamics,  elasticity,  oil  reservoir  simulation  and 

other  areas  of  engineering  when  diffusion  is  neglected.  In  solutions  of  (1), 

information  travels  at  speeds  and  X2  given  by  the  eigenvalues  of  the 

3f 

matrix  — .  Since  this  matrix  depends  on  u,  v,  the  speeds  X^  and  X2 

depend  on  the  solution,  and  this  leads  to  the  formation  of  discontinuities 

called  shock  waves.  The  present  paper  deals  with  the  classification  of 

Rlemann  problem  solutions  (solutions  which  evolve  from  a  single  discontinuity 

3f 

A# 

at  tljae  t  ■  0)  near  an  isolated  point  at  which  t—  is  a  multiple  of  the 
identity,  so  that  X^  >  X2.  Such  a  singularity  has  no  analogue  in  the  linear 
theory. 
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THE  CLASSIFICATION  OF  SOLUTIONS  OF  QUADRATIC  RIEMANN  PROBLEMS  (I) 


E.  Isaacson**'*'^'®' D.  Marchesin***^'^'^,  B.  Plohr^'®*’*®  and  B.  Temple’ ' ’'® 


§1 .  Introduction 

We  consider  the  general  2x2  system  of  conservation  laws  with  quadratic  flux 
functions 


(1) 


Uj.  +  {a^u^  +  2b^uv  +  -  0  , 

’'t  *  1  *  ^*’2“''  * 


0  . 


'2’  ^x 

In  particular,  we  solve  the  Rlemann  problem  globally  for  a  specific  range  of  the 
coefficients  a^,  b^^,  c^- 

System  (1)  is  of  interest  because  solutions  of  (1)  approximate  solutions  of  an 
arbitrary  2x2  system  of  conservation  laws 


(2) 


Ht  *  £*~’x  ~  ' 


where  u  =  (u,v),  f  =  (f,g),  in  a  neighborhood  of  an  Isolated  hyperbolic  singularity. 
Such  a  singularity  is  an  isolated  point  in  a  neighborhood  of  which  (2)  is  hyperbolic  and 
at  which  the  Jacobian 


A(ii) 
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has  equal  eigenvalues  and  Is  dlagonalisable.  In  fact,  systea  (1)  Is  obtained  froa  systea 
(21  as  follows-  Let  <  k^(a)  denote  the  eigenvalues  of  A(u),  and  let  u^  denote 

the  Isolated  point  at  which  ^  ^0'  replace  u  by  “  "  “q  •“<* 

translate  the  reference  frame  (x,t)  to  (x  -  Xgt.t)  so  that  the  resulting  systea  has  an 
Isolated  singularity  at  u  •  (0,0)  with  corresponding  double  eigenvalue  \  •  0.  Then 
system  (1)  Is  obtained  by  expanding  the  flux  functions  of  this  transformed  system  in 
Taylor  series  about  (0,0)  and  neglecting  higher  order  terms. 

System  ( 1 )  can  be  reduced  further  by  a  nonsingular  linear  change  of  dependent 
variables.  Two  systems  related  by  such  a  transformation  S  are  isoiK>rphlc  in  the  sense 
that  Su(x,t)  Is  a  weak  solution  of  the  transformed  system  If  and  only  If  u(x,t)  Is  a 
weak  solution  of  the  original  system.  Since  the  nonsingular  transformation  S  contains 
four  free  parameters  and  since  system  (1)  contains  six  parameters,  we  expect  to  find  a  two 
parameter  family  of  Isomorphism  classes  for  system  (1).  Thus,  we  look  for  representatives 
of  the  Isomorphism  classes  In  a  normal  form  containing  two  free  parameters  [8] .  In  what 
we  consider  to  be  a  breakthrough.  Shearer  and  Schaeffer  showed  In  [17]  that  when  system 
(1)  Is  hyperbolic,  there  Is  a  nonsingular  linear  change  of  dependent  variables  which 
transforms  system  ( 1 )  Into 

1  5  2 

u^  —  {au*  +  2buv  V  »  0  , 

(3)  , 

v^  +  —  {bu  +  2uv}jj  ”  0  . 

System  (3)  depends  on  two  free  parameters  a  and  b  and  can  be  taken  as  a  normal  form 
for  the  hyperbolic  quadratic  systems  (1).  It  is  also  shown  In  [3,17]  that  the  Integral 
curves  of  (3)  fall  Into  four  nonisomorphic  classes  depending  on  the  parameters  a  and  b. 
These  classes  define  four  regions  In  the  a,b-plane  which  are  referred  to  as  Regions  1-IV. 

The  regions  are  determined  by  the  number  of  lines  which  form  the  Hugonlot  locus  of 
the  origin,  as  well  as  the  direction  of  Increase  of  the  appropriate  eigenvalue  on  these 
lines.  In  Regions  I-III,  the  Hugonlot  locus  consists  of  three  distinct  lines,  while  in 

Region  IV  it  consists  of  one  line.  Specifically,  the  boundary  between  Regions  I  and  II  Is 

3  p  2 

given  by  a  *  —  the  boundary  between  Regions  II  and  III  is  given  by  a  «  1  +  b  #  and 


.«  *-•  •jK 


the  boundary  between  Regions  III  and  IV  Is  given  by  4{4b^  -  3(a  -  2)}^  -  {16b^ 

*  9(1  -  2a)b}^  (cf>  17].  The  structure  of  solutions  in  each  region  is  simplest  when 
b  •  0  since  then  solutions  have  both  up-down  synmetry  [(u,-v)  satisfies  (3)  if  and  only 
if  (u,v)  does]  and  left-right  syaanetry  ( (u(x,t),v(x,t) )  satisfies  (3)  if  and  only  if 
(-u(-X(t),v(-x.t) )  does].  We  call  the  systems  with  b  •  0  symmetric.  (An  additional 
siiqplifying  feature  of  the  symmetric  systems  is  that  the  lines  on  which  genuine 
nonlinearity  falls  coincide  with  the  Hugoniot  locus  of  the  origin.) 

The  present  paper  is  the  first  of  a  series  in  which  we  give  the  solution  of  the 
Riemann  problem 


X  <  0 

(4) 

u(x,0  )  »  j 

1  "r  =  ‘“r'V' 

x>  0 

for  the  symmetric  systems  In  Regions  II-IV  of  Schaeffer  and  Shearer  [17].  There  are  new 
features  in  these  regions  that  do  not  occur  In  Region  I  [cf.  18].  First,  a  type  of  shoc)c 
which  we  call  coBg>rea8ive  appears  in  solutions  of  the  Riemann  problem,  and  the  existence 
of  such  shocks  is  necessary  to  ensure  (in  the  x,t-plane)  the  continuous  dependence  of 
solutions  on  u  and  u  [cf.  6,7,9,21].  Schaeffer  [25]  noted  that  the  compressive 
shocks  perturb  to  a  1-shock  followed  by  a  2-shock.  This  fact  is  manifested  in  the  triple 
shocks  which  appear  in  the  solutions. 

A  second  feature  that  occurs  in  Regions  II-IV  that  does  not  occur  in  Region  I  is  the 
following:  in  Regions  II-IV  there  are  lines  on  which  X  “  0,  whereas  in  Region  I,  X  •  0 

if  and  only  if  u  ■  0.  These  lines  play  a  central  role  both  in  the  classlf Ication  as  well 
as  in  the  structure  of  Riemann  problem  solutions.  In  fact,  the  lines  \  =  0  give  a 
geometric  interpretation  for  the  boundaries  between  the  Regions  I-III  located  by  Schaeffer 
and  Shearer.  (See  the  Appendix  and  Figure  6.)  Also,  in  Regions  II-IV  the  structure  of 
the  shock  types  on  the  Hugoniot  locus  of  a  point  u  changes  as  u  crosses  a  line  X  •  0; 
this  entails  a  corresponding  change  in  the  structure  of  the  solutions  of  the  Riemann 
problem.  The  reason  for  the  change  in  shock  type  is  that  the  state  u  on  the  line  X  *  0 
can  be  joined  to  the  state  -u  by  a  two-sided  contact  discontinuity. 
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Our  construction  of  the  solutions  of  the  Rieaann  problen  is  based  on  a  nuaerlcal 


construction  of  the  Hugoniot  loci  together  with  the  structure  of  the  integral  curves 
obtained  from  [171.  Some  analysis  is  presented  to  justify  the  salient  features  of  the 
solutions  obtained.  In  particular,  our  analysis  uses  an  explicit  paraMterlsation  of  the 
Hugoniot  locus  (see  the  Appendix).  The  fact  that  the  Hugoniot  locus  is  star-like  with 
respect  to  the  left  state  was  pointed  out  by  Shearer  [24]. 

In  the  present  paper  we  give  the  solution  of  the  Riemann  problem  for  system  (3)  in 
the  parameter  range 


(5)  a  >  2,  b  •=  0  . 

This  condition  specifies  the  symmetric  systems  in  Region  IV.  We  present  the  solution  in 
Section  3  by  means  of  a  series  of  diagrams  (Figures  7A-F).  Because  of  up-dom  symmetry, 
we  give  the  solutions  only  for  u^  in  the  lower  half-plane.  The  solution  diagrams  are 
qualitatively  the  same  for  all  u^  in  the  sector 

A,  =  {u  :  6^  <  6^  <  0)  , 
and  are  qualitatively  the  same  for  all  u  in  the  sector 

A2  »  {u  :  -»  <  <  e,)  , 


where 


9.  *  arctanf — I 
L  Vu^^ 


0^  =  arctan(-/a)  . 

This  separation  angle  8*  depends  on  a,  and  the  ray  6  ■  is  the  ray  (in  the  lower 
half-plane)  on  which  =  0.  The  solution  consists  of  a  1-composite  wave  followed  by  a 
2-composite  wave  in  analogy  with  the  local  solutions  for  strictly  hyperbolic  systems  [cf. 
11,14].  Here,  however,  the  wave  curves  have  a  different  structure:  as  in  [6,7,9,21], 
intermediate  solution  states  do  not  depend  continuously  on  the  data  due  to  the  appearance 
of  compressive  shocks.  However,  as  in  (6,7,9,21),  continuous  dependence  on  u  and  u 
is  ensured  in  x,t-space  because  of  the  coincidence  of  shock  speeds  in  the  compressive 


shocks.  In  this  paper  we  construct  explicitly  for  each  state  u  the  1-wave  curve  ^^(u) 
and  the  2-wave  curve  1^2 such  that  the  following  thsoraa  holds: 

THEOREM:  For  each  pair  of  states  u.  and  u_.  there  exists  an  Interaedlate  state 
-M  ^  such  that  U|^  ;  solution  of  the  Riemann  problem  (3),  (4), 

(5)  consists  of  the  lixave  froai  Uj^  to  followed  by  the  2-wave  from  u^  to 

Moreover,  the  solution  Is  unique  In  x,t-space  and  depends  continuously  on  u^  and 


In  Figures  7a-F,  pictorial  solutions  of  the  Riemann  problem  are  presented  In  which 
the  state  u^  Is  fixed,  and  an  arbitrary  point  in  the  diagram  represents  u^>  The  waves  in 
the  solution  of  the  Riemann  problem  are  determined  by  the  path  from  to  u^^  which 

consists  of  (portions  of)  wave  curves  and  is  indicated  by  the  arrows. 

In  Section  2  we  discuss  the  wave  curves  (tff  and  ^2’  Section  3  we  describe  the 
diagrams  in  Figures  7A-E  individually,  ^-'.d  in  Section  4  we  verify  the  salient  features  of 
these  diagrams.  General  properties  of  system  (1)  are  stated  in  the  Appendix  and  are 
referred  to  throughout.  The  Appendix  also  contains  the  interpretation  of  Regions  I-III  In 


terras  of  the  lines  A  «  0.  These  are  depicted  in  Figures  8A-D 


The  Integral  curves  of  the  eigenvector  fields  of  the  Jacobian  A()})  detemlne  the 
rarefaction  waves  of  system  (2).  The  Hugonlot  loci  H(a^)  for  states  determine  the 

shock  waves  of  system  (2),  and  are  given  by  [20] 

M(Uj^)  =  {u  :  s(u  -  Uj^)  «  ftjj)  "  some  s  c  R  }  . 

The  general  solution  of  the  Rlemann  problem  which  we  construct  is  obtained  by  coiq>oslng 
rarefaction  waves  and  shock  waves.  Because  the  flux  functions  In  (1)  are  homogeneous 
quadratics,  the  Integral  curves  and  Hugonlot  loci  for  system  (1)  have  the  following 
scaling  property:  if  T  is  an  Integral  curve  or  Hugonlot  locus  through  u, ,  then  cF 

'“la 

is  the  corresponding  Integral  curve  or  Hugonlot  locus  through  cu^^,  c  0.  The 
eigenvalues  and  shock  speeds  scale  similarly.  Thus  the  integral  curves  and  Hugonlot  locus 
through  u^  determine  the  corresponding  Integral  curves  and  Hugonlot  locus  for  each  state 
u  satisfying  9  =  6^^  where 

B  “  arctan('^)  . 

For  the  Hugonlot  loci,  the  up-down  and  left-right  symmetries  Imply  that  If  u^^  c  (-/(Uj^), 
then  u„  f  H(u  )  and  u  r  f((u  )  where  u  -  (u,-v)  and  u  ■  (-u,v/.  'Rius,  the 
reflection  of  a  Hugonlot  curve  about  either  the  u-  or  v-axis  Is  also  a  Hug^nUt 
curve.  In  terms  of  wave  speeds  and  shock  speeds,  the  symmetries  of  (3)  when  b  ■<  0  yield 
the  following  relationships: 

(2.1)  X  (u)  =•  X  (u)  -  -X-(ii),  p  =  1,2 

(2.2)  “  ^^~L'~R* 
where  p  denotes  the  other  family. 

In  Table  1  we  label  nine  types  of  shocks  according  to  the  Inequalities  that  hold 
between  the  characteristic  speeds  of  u^^,  Uj^  and  the  shock  speed  ci(Uj^,Uj^).  (Obvious 
inequalities  are  omitted.)  Since  the  two  wave  speeds  X^  and  X2  assigned  to  each  state 
cjre  ordered.  Table  1  gives  the  totality  of  shock  types  that  can  occur  in  any  2  »  2 
system  (2).  In  the  case  of  the  symmetric  systems  (3),  the  formulas  (2.1)  and  (2.2)  imply 
the  following  symmetries  for  shock  types:  If  Uj^  e  H(Uj^),  then  the  shock  type  of 


•r.  •r.  .  .  .  .  .  . 


<Uj^,Uj^>  will  be  the  same  as  the  shock  type  of  <Uj^,Uj^>»  and  the  shock  type  of 
will  be  what  we  call  the  inverse  of  the  shock  type  of  <u, /U_>.  Here,  the  inverse  shock 
type  is  obtained  by  reversing  the  inequalities  which  hold  between  the  shock  speed  and  each 
wave  speed  to  the  right  and  left  of  the  shock  and  also  by  interchanging  the  families  1  and 
2.  Therefore  [cf.  Table  1],  1-shock  and  2-expanslve  are  inverses;  2-shock  and  1-expansive 
are  inverses;  compressive  and  expansive  are  inverses;  right  transport  and  left  transport 
are  inverses;  and  crossing  is  its  own  inverse.  Thus  in  the  symmetric  cases,  all  Hugoniot 
curves  and  shock  types  are  determined  by  the  Rugoniot  curves  and  shock  types  for  states  in 
a  single  quadrant  of  the  u,v-plane.  Among  the  above  shock  types,  some  occur  in  the 
solutions  of  system  (3)  that  do  not  occur  for  strictly  hyperbolic  systems.  For  the 
symmetric  systems  of  classes  II-IV,  there  are  three  types  of  shocks  that  appear  in 
solutions:  l-shocks,  2-shocks,  and  compressive  shocks  which  are  defined  by  [cf.  Table  1] 

We  refer  to  any  of  the  three  types  of  shocks  above  as  admissible.  We  refer  to  shocks  of 
other  types  as  inadmissible.  (We  note  that  the  transport  shocks  are  the  only  types  that 
do  not  appear  in  Hugoniot  loci  of  the  symmetric  systems  in  classes  II-IV.)  Our  scheme  for 
labeling  shock  types  in  Figures  3-7  is  described  in  the  legend  which  precedes  the 
diagrams. 
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*-n  the  didgrams  above,  the  centerline  denotes  a  shock  in  x,t-space.  The  lines  to  the  left 
and  right  of  the  shock  are  characteristic  lines,  and  indicate  the  relationship  between  the 
shock  speed  and  the  wave  speeds  on  the  left  and  right. 


f2.1.  Integral  Carve 

For  ayatea  (3)  with  paraaetar  valuaa  (S), 


’aa  V 

/ 

_  V  u_ 


and  the  Integral  curvea  are  depicted  In  Figure  1  [3,17].  He  call  the  Integral  curvea 
associated  with  the  eigenvalues  lf(u)  <  2-lntegral  curves.  For  exaaple, 

when  a  »  3  the  Integral  curves  are  the  one-paraaeter  faally  of  parabolas 

where  Ug  Is  the  value  at  which  the  parabola  crosses  the  u-axla.  The  u-axls  is  itself  an 
Integral  curve  for  the  aysHMtrlc  aysteas  b  *  0.  For  general  a  >  2,  the  Integral  curves 
are  parabola-llke  and  the  t-lntegral  carves  open  to  the  left,  the  2-lntegral  curves  open 
to  the  right.  The  arrows  in  Figure  1  Indicate  the  direction  In  lAkich  the  corresponding 
eigenvalues  Increase.  He  define  the  l>rarefection  curve  1  ~  <>2,  to  be  the  set 

of  states  ^  such  that  the  solution  of  the  Rleaann  problaa  (3),  (4)  is  a  pure  1- 
rarefsctlon  wave  (see  (14});  connected  portion  of  the  l-lntegral 

curve  through  u^  consisting  of  those  states  g  for  which  Increases  along  the  curve 
froa  ^  to  u. 


{2.2.  Ihe  Buoonlot  Locus 

For  systea  (3)  with  paraaeter  values  (5),  the  qualitative  shape  of  the  Rugonlot  locus 
H(Uj^)  Is  as  follows:  for  Uj^  •  S»  W(Uj^)  Is  the  u-axls;  for  Uj^  on  the  u-axls,  H(Uj^) 
consists  of  an  ellipse  together  with  the  u-axls;  and  for  u^  off  the  u-axls,  H (u^) 
perturbs  to  a  closed  loop  surrounding  the  origin  with  two  tails  which  are  asyiqitotlc  at 

Infinity  to  opposite  ends  of  the  horlsontal  line 

a 

V  - - -  V.  . 

a  -  2  !• 

These  properties  follow  froa  the  explicit  paraaeterltatlon  of  the  Hugonlot  locus  given  In 
the  Appendix  and  are  depicted  In  Figures  3A-F.  The  shock  type  Is  also  indicated.  He 


let  S^(a^)  denote  that  portion  of  coneleting  of  1-shocks  for  p  •  If  2-shocks 

for  p  »  2,  and  ccnpresslve  shocks  for  p  •  c. 

The  subset  Is  non-empty  for  In  the  sector  A^,  and  Is  empty 

for  In  the  sector  A2>  Moreover f  S^(tt^)  has  two  disconnected  components  for 

e  A^  but  only  one  for  u^  e  ^2’  Ais  explains  «ihy  the  ray  6-8*  is  the  boundary 
across  which  the  qualitative  structure  of  the  solution  dlagrasw  changes* 

We  note  in  Figure  3B  that  Bj^,  are  the  points  at  which 

in  Figures  3A,  B,  C,  D  that  Cj^  is  the  point  at  which 

and  in  Figure  3C  that 


S2.3.  Composite  Waves 

The  solutions  of  the  Riemann  proble;|  are  cospositions  of  pure  rarefaction  waves  and 
pure  shock  waves.  We  refer  to  a  Riemann  prob}.em  solution  consisting  entirely  of  waves  of 
one  family  as  a  composite  wave  of  that  family.  The  composite  waves  in  the  solutions  of 
Riemann  problems  for  system  (3)  in  parameter  range  (5)  are  of  two  types:  a  1 -rarefaction 
wave  followed  by  a  1 -shock  wave  and  a  2-shock  wave  followed  by  a  2-raref action  wave.  In 
each  case,  the  speed  of  the  shock  equals  the  speed  of  the  adjacent  rarefaction  wave.  We 
call  these  IRS-  and  2SR-coaposite  waves,  respectively.  (See  Figure  4.  nie  letters  in 
Figure  4  refer  to  states  appearing  in  Figure  5  where  the  composite  waves  may  be  viewed  in 
state  space.) 

For  each  u  we  define  the  1-composite  curve  C^(u).  In  Figure  5  a  IRS-composite 
wave  with  left  state  u  is  a  wave  consisting  of  a  1-rarefaction  wave  from  u  to 
F  (  R^(u)  followed  by  a  1-shock  from  F  to  D  e  S^(F)  such  that  a(F,D)  “  X^(F).  Such 
a  shock  exists  if  and  only  if  F  lies  on  that  portion  of  ^^(v)  ii'  Sector  A^.  We 


define  the  1-composite  curve  C, (u)  to  be  the  set  of  all  such  states  D.  Thus,  the  1- 


coiiposite  curve  C^(u«)  ia  aaxlaal  for  on  the  boundary  ray  6  <•  6*  in  the  sense 

that 


-  C,(u) 

for 

u«  e  R^(u)  , 

Ci<S.) 

2 

for 

u  e  R^(u,)  ; 

for  exaaplSf  in  Figure  5, 


C^(u)  -  C,(u,)  -  tEDC.l  and  C,(F)  -  [ED)  C  C,(u«)  . 

(Here,  letters  enclosed  by  brackets  denote  the  curves  that  connect  the  corresponding 
points  In  the  Figures*)  In  fact,  as  F  noves  fron  E  to  u^  along  R^(u,,),  0  moves 
from  E  to  C*  along  C^(j2,  )•  Moreover, 


C«  “  -u*  , 


and 


-  XjCC*)  -  X,()J,)  “  0  . 

For  each  curve  T  we  define  the  3-boundary  curve  based  on  r*  In  Figure  5,  a  2SR- 

coi^>oslte  wave  with  left  state  u  Is  a  wave  consisting  of  a  2-Bhock  wave  from  u  to 

C  «  ^(u)  followed  by  a  2-rarefaction  wave  from  C  to  any  G  c  such  that 

<r(u,C)  •  a  shock  exists  if  and  only  If  u  lies  on  that  portion  of  r  In 

sector  ^2’  He  define  the  2-boundary  curve  to  be  the  set  of  all  such  states  C.  For 

example,  in  Figure  5,  [CeC*]  Is  the  2-boundary  curve  based  on  (u,uw].  Thus,  for  left 

states  u  on  T  »  {“.a"!*  the  states  C  on  the  2-boundary  curve  (C*C»]  are  transition 
states  between  2-shock  waves  and  2-coiig>osite  waves* 


j2.4.  Wave  Curves 

He  define  the  wave  curve  (</^(u^),  1  ~  1,2,  to  be  the  set  of  states  Uj^  for  which 
the  solution  of  the  Riemann  problem  (2),  (4)  Is  a  pure  composite  wave  of  the  appropriate 
family*  nius,  <see  Figures  6A,  B) 

and 
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i3.  Solution  of  the  Rlemann  Problem 

For  each  u^  we  describe  the  solution  of  the  Rlemann  problem  for  arbitrary  u^  by 
means  of  a  diagram.  In  Figures  7B  and  D  we  present  a  diagram  of  the  solution  of  the 
Rlemann  problem  for  representative  values  of  u^  In  sectors  and  Aj  In  each  of 

which  the  solution  diagrams  are  qualitatively  the  same.  For  completeness  and  in  order  to 
see  clearly  the  continuous  dependence  of  solutions  on  u, ,  we  also  give  in  Figures  7A,C,E 

'“la 

the  solution  diagrams  for  typical  values  of  on  the  boundaries  of  sectors  A^  and 

Aj.  i.e.,  for  Uj^  on  the  rays  0  «  0,  0  •  0*,  and  B  ”  -v.  In  Figures  7a-D,  the 
representative  values  u  are  chosen  to  lie  on  the  1-lntegral  curve  through  the  fixed 

state  E  on  the  positive  u-axls.  In  Figure  7F  we  present  the  solution  for  u  «  0. 

The  solution  of  the  Rlemann  problem  consists  of  a  1-wave  with  left  state  u  and 

'“la 

right  state  u^^  followed  by  a  2-wavc  with  left  state  u^  and  right  state  u^^.  We  find 

the  Intermediate  state  In  the  Figures  as  follotvs:  given  follow  the  2-wave  curve 

backwards  from  u.  (opposite  the  direction  of  the  arrows)  until  you  reach  a  point  u  in 
the  1-wave  curve  W,(u, ).  The  state  u  so  constructed  satisfies  u  c  W  (u  )  and 
defines  the  waves  in  the  solution.  This  procedure  Is  not  well-defined  when 
form  a  triple  shock; 

(3.1)  u,^  f  Sr  *  “r  *  ' 

and 

(3.2)  "  °'!iM'HR'  • 

Triple  shocks  occur  when  u_  f  S  (u  )  or  when  u  Is  In  the  triple  shock  curve  (which 

*“R  C  '^la  '“R 

will  be  defined  In  each  Figure).  When  the  procedure  Is  not  well-defined,  there  are  two 
such  Intermediate  states  u  ;  however,  the  two  solutions  obtained  are  identical  in  x,t- 
space  because  then  all  shock  speeds  In  the  problem  are  equal.  This  ensures  continuous 
dependence  of  the  solution  on  the  Initial  data  and  u^^. 

We  now  discuss  the  solution  diagrams  7a-E  Individually. 
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j3.1.  Solution  for  6^  “  0 

In  Figure  7a,  i  E  lies  on  the  y-axls  to  the  right  of  the  origin.  In  this  case 
the  Hugonlot  locus  consists  of  the  jj-axls  together  with  the  ellipse  E 

A 

depicted.  The  vertices  of  E  on  the  major  axis  are  u^^  and  ^  Bie 

point  A  is  given  by 


Point  A  is  the  limit  of  intersections  of  S2(u)  with  the  u-axls  as  u  tends  to 
through  states  with  v  ^  0.  The  solution  of  the  Riemann  problem  consists  of  a  1-shock 
with  left  state  u^  and  right  state  u^^  £  S^{Uj^)  followed  by  a  2-rarefaction  wave  if 
u  lies  outside  E,  and  followed  by  a  2-shock  if  u-  lies  inside  E.  n>e  u-axls 
between  Cj^  and  A  is  S^(u^)  where  triple  shocks  occur  in  pairs.  In  fact,  for  each 
V  £  E,  V  /  Uj^,  there  exists  a  unique  jj  e  such  that  both  {Uj_,v,u)  and 

fu  ,v,u)  are  triple  shocks  of  the  same  speed.  (Here  v  is  the  reflection  of  v  in  the 
u-axis.)  Thus  for  u^^  c  S^(Uj^),  the  solution  of  the  Riemann  problem  is  unambiguous  in 
the  x,t-pl«me. 

§3.2.  Solution  for  0»  <  <  0 

In  Figure  7b,  Uj^  lies  in  the  sector  A.|.  The  Hugonlot  locus  depicted  by 

the  thin  solid  curve  through  u^^.  Here 

Si'sEl'  -  ^  ' 

C,(Uj^)  -  (EDj^l  , 

and  the  union  of  these  three  curves  is  the  1-wave  curve  (;/^(u^).  Moreover, 

and  the  triple  shock  curve  is  (3^*^ ' 

For  Ujj  on  either  of  these  latter  two  curves,  triple  shocks  occurs  for 

R  f  ^  (u, )  there  is  a  pair  of  triple  shocks  {u, ,P,R)  and  {u, ,Q,R}  with  all  shock 
c  '^L  « 


speeds  belt!?  equal*  for  R  e  [B^A] ,  a  single  triple  shock  occurs.  Hence,  continuous 
dependence  Is  ensured. 

The  curve  [Cj^C'“J  Is  the  2-boundary  curve  based  on  [KjU'*],  and  2SR-coaiposlte 
waves  occur  precisely  when  Uj^  lies  "above"  [C^C**]  so  that  e  [KjU'“]  . 

The  states  and  0^  satisfy 

The  state  Is  the  point  where  H(Uj^)  la  tangent  to  a  2-lntagral  curve.  (See  the 

Tangency  Rule  In  the  Appendix.)  This  Implies  that  a(Uj^,Cj^)  «  that  Is 

In,  and  marks  the  end  of,  the  2-boundary  curve  (Cj_C‘«>].  This  Insures  continuous 
dependence  near  the  point  C^. 

For  In  sector  A^,  the  states  B^,  separate  the  states  Uj^, 

portion  of  consisting  of  the  closed  loop.  As  decreases  to  6*,  Bj_  and 

become  coincident  with  Cj^,  and  so  both  *****  portion 

of  disappear.  The  fact  that  there  exists  a  point  at  which  H(Uj_)  Intersects  the 

l-composlte  curve  characterizes  the  states  In  A^. 

^3.3.  Solution  for  “  9* 

In  Figure  7C,  u^^  “  u,  lies  on  the  line  0  -  0*.  The  Hugonlot  locus  H(u^)  Is 

depicted  by  the  thin  solid  curve  through  u^^.  Here 

S^(Uj^)  -  ' 

and 

The  union  of  these  three  curves  Is  U/^(Uj^).  In  this  case 

' 

and  the  curve  [C,Al  Is  the  triple  shock  curve. 

The  curve  lC*C'"I  is  the  2-boundary  curve  based  on  (Uj^'*),  and  2SR-composlte 
waves  occur  precisely  when  lies  "above"  (C*C*“]  so  that  u^^  £  *  Note  that 


^ 'A '5'. A 


C*  >=  -u  is  the  point  where  H(u  )  is  tangent  to  a  2-rarefaction  curve,  and  also  that 

a(^,C,)  »  X^(Uj^)  -  ^^(C,)  . 

Thus  9  •=  8,  is  the  ray  where  B^,  Cj^  and  are  coincident. 

§3.4.  Solution  for  -s  <  0^  <  9^ 

In  Figure  7D,  u^  lies  in  sector  A2.  the  Hugoniot  locus  is  depicted  by  the 

thin  solid  curve  through  u^.  Here 

and 

C^(u^)  =  tEDC*l  . 

The  union  of  these  three  curves  is  Xn  this  case 

(3.4)  ^  ' 

and  [ABC*]  is  the  triple  shock  curve.  Condition  (3.4)  characterizes  the  states  in  A2. 

For  U|^  e  [u.EC*],  the  solution  consists  of  a  1-rarefactlon  wave  ta)clng  u^^  to  u* 
followed  by  the  solution  of  the  Rlenann  problem  <u,,Ujj>  given  in  Section  3.3.  The 
curve  (C«Cj_«»]  is  the  2-boundary  curve  with  base  curve  again  2SR-coaposite 

waves  occur  precisely  when  u  lies  "above"  the  curve  (C*C,«)  so  that  u  e  ("U  u.J. 

— R  ^(4  .*!/*• 

Note  that  C*  Is  the  end  of  the  1 -composite  curve  in  the  sense  that  there  do  not 
exist  1-shocks  with  left  states  u  f  which  have  speed  X^(u). 

§3.5.  Solution  for  9j^  =  -it 

In  Figure  7E,  u^^  lies  on  the  negative  u-axls.  The  Hugoniot  locus  consists  of  the 

u-axis  together  with  the  ellipse  E  depicted.  The  vertices  of  E  on  the  major  axis  are 
u  and  G  =  -( — — r]  u  ,  and 

L»  3  ~  £  '*'Li 

and 


-16- 


§4.  Justif icatlon 

The  solution  of  the  Riemann  problem  presented  in  Section  3  is  based  on  the 
assumptions  that  the  integral  curves  and  Hugoniot  loci  have  the  qualitative  features 
discussed  in  Section  2.  The  relevant  features  of  the  Integral  curves  which  are  required 
for  our  analysis  are  verified  in  13,17],  The  features  of  H(u-)  (including  shapes  and 
locations  of  shock  types)  have  been  verified  numerically,  using  an  explicit 
parameterization  of  (see  the  Appendix).  We  here  present  analytical  evidence  which 

supports  the  results  of  these  numerical  computations.  Note  that  for  u  off  the  u-axls, 
system  (3)  in  parameter  range  (5)  is  strictly  hyperbolic  and  genuinely  nonlinear  in  a 
neighboihood  of  u^^,  and  so  the  local  theory  of  Lax  (11)  applies  in  some  neighborhood  of 
.  o'e  now  discuss  the  global  features  of  the  Hugoniot  loci  depicted  in  Figures  3A-C. 


Figure  3A:  In  this  case  u^  is  on  the  positive  u-axis.  The  Hugoniot  locus 
consists  of  the  u-axis  together  with  the  ellipse  E  given  by 


with  renter 


(a  -  2)  ~L 


(a  -  2)u  +  2UlU  +  V  =  au£  , 

For  u  =  (u,0),  we  have  (see  the  Appendix) 


1 


o(u,,u)  =  —  a(u  +  Ut)  , 
--L  ~  2  L 


( 1  ^  (  u)  ,  X^lu) }  =  {au,u}  and  {X  ^  (u^^ ) ,  X ^  *  (*“L'“l^  ' 
whicl.  verifies  the  shock  types  on  the  axis,  and  verifies  that  the  shoc)t  type  changes  at 


a  ...  a  -  2 

u  0  ana  As"*  ■■ 


a  -  2  -“-L 


Figure  3B :  In  this  case  f  A  .j .  The  points  and  are  points  where  a 

transition  in  shock  type  occurs.  The  existence  of  the  point  where  tangent 

to  a  2-integral  curve  follows  from  the  geometry  of  the  integral  curves.  This  tangency 
dt  implies  (see  the  Tangency  Rule  in  the  Appendix) 


Nf)te  also  that  u_  and  -u_  both  lie  on  H(u_)a  and 


This  1*  obtained  from  the  Midpoint  Rule  or  directly  froa  the  Ranklne-Hugonlot  condition 


(see  the  Appendix).  Also  the  speed  X^(u^)  satisfies 


°  < 


a  +  1 


~  /r -  t  4v? 


<  “l  ■  ■’‘"l'Hl'  ' 


2  “L  2  ■  ■'  T.  "  L 

where  the  left  hand  Inequality  is  equivalent  to  u  e  A^.  Ilius  the  shock  speed  a(u^/U^) 

Is  larger  than  X^(u^),  but  o(Uj^,-u^)  la  snvaller  than  X^(Uj^);  hence,  there  exists  a 


point  r  W(u^)  between  u^  and  at  which  o(Uj^,D^)  »  X^{Uj^).  Similarly, 


Is  the  Intersection  of  H(u  )  with  the  negative  u-axls,  and  satisfies  (see  Appendix) 

***Li 

Therefore,  there  Is  a  point  above  the  u-axls  between  -u^^  and  u^  at  which 

o(Uj^,Bj^)  •  Xj(u^).  Consequently,  *  triple  shock. 


Figure  3C;  In  this  case  -  ji,  la  on  the  ray  8-8*.  As  decreases  to  8,  in 

Figure  3B,  S^(aj^)  the  portion  vanish.  In  fact, 

indicating  that  the  points  B^  and  do  Indeed  become  coincident  with  -u,  as  9j^ 

decreases  to  0«-  We  also  note  that  X2(-u*)  -  o(u,,-u,)  implies,  by  the  Tangenoy  Rule 
that  W(u,)  is  tangent  to  a  2-integral  curve  at  -u«  and  thus  C,  “  -u, . 


APPENDIX 


Here  we  list  for  reference  several  properties  of  the  quadratic  conservation  laws  (1> 
which  are  helpful  in  the  study  of  the  Riemann  problen- 


MIDPOINT  RULE;  For  system  (1),  Hp  t  line  segment  joining  u^^ 

to  Uj^  is  tangent  to  a  p-integral  curve  at  the  midpoint  of  the  segment.  Moreover,  in 


this  case 
(A.1) 


3f 

Proof:  Use  the  linearity  of  —  together  with  the  general  formula 

■  dU 


1  3f 


s(u  -  u  )  =  f(u  )  -  f{u  )  =  {/  T—  (U  ♦  t(u  -  u  ))dT)*(u  -  u  ) 

~R  ~L  ~  ~R  ~  ~L  ^  3u  ~L  ~R  ~L  ~R  ~L 


u^  f  '  S2  *  '  ^3  ^  • 


TRIPLE  SHOCK  RULE;  For  system  (2),  suppose  the  states  u^,  u^,  and  u^  satisfy 
{A. 2) 

Then  either  u.,  u,  and  u,  are  colinear  or  else 
~2  ~'3 

(A. 3)  o(u2,u^)  •  oCUj.Uj)  •  a(u^,Uj)  . 

Proof;  Let  s^j  =  a(Uj^,u.).  Then 


=2l'!i2 

-  s,) 

»  flUj) 

■  ' 

®32<!i3 

■  S2> 

-  f(Uj) 

-  f{U2)  , 

®13'Sl 

-  a3> 

=  f(u,) 

-  ftUj)  . 

Adding  gives 

(S21 

8,3>(U2  -  u,)  t 

‘®32 

-  s,3)(u 

3  -  “2’  '  °  ' 

and  the  result  follows. 

TANGENCY  RULE:  For  system 

(3),  assume  that 

5iL 

1  H(0). 

Then  the  following 

equivalent  regarding  u  t  H(u^),  u  ^  u^^: 

(i)  tangent  to  a  p-integral  curve  at  the  point  u. 


(11)  8  ~  0  at  u.(Dot  denotes  differentiation  with  respect  to  arc  length  along  H(u-).) 

( ill )  X  (u)  “  a(u, ,u) . 

P  ~  ~L  ~ 


Proof:  Differentiate 


along  W(u^)  to  obtain 


i(u  -  Uj^)  ”  f(u)  - 


8(u  -  U-)  +  au  “  A(u)u  , 


[A(u)  -  sl]u  “  s(u  -  u  )  . 


Suppose  H(^)  is  tangent  to  a  p-lntegral  curve.  Then  u  is  a  right  eigenvector  of  A 


with  eigenvalue  X^  so  that  (A. 4)  becomes 

(A. 5)  (X  -  8)a  “  s(u  -  u.)  . 

p  ~  ~  ~L 

Since  Uj^  ^  line  joining  u^^  to  u  is  not  tangent  to  H<Uj^)  at  u.  Thus  we 

conclude  from  (A. 5)  that  s  “  0  and  s  ”  Xp,  so  that  (1)  “>  (11)  and  (1)  ->  (ill).  Now 
suppose  (11).  Then  (A. 4)  gives  that  u  Is  en  eigenvector  and  s  «  Xp  for  p  “  1  or  2. 
Tbus  (11)  «>  (1)  and  (11)  “>  (ill).  Finally  suppose  (ill):  s  •  Xp(u).  If  u  Is  not 
parallel  to  the  eigenvector  of  Afji)  and  i  i*  0  then,  (A. 4)  Implies  that  u  - 

Is  parallel  to  r^,  (the  other  eigenvector  of  A(^)).  However,  by  the  midpoint  rule, 
u  -  Uj^  Is  also  tangent  to  an  integral  curve  at  -j  (Uj^  +  u)  .  But  this  contradicts  the 
fact  that  If  Uj^  ^  then  every  line  through  is  tangent  exactly  once  to  an 

Integral  curve.  Consequently,  u  is  tangent  to  ^  at  u,  and  hence  s  “  0.  Thus 
(ill)  “>  (1)  and  (ill)  ->  (11). 


EIGENVALUES  AND  SHOCK  SPEEDS:  For  system  (3)  with  b  -  0,  the  eigenvalues  are 

X^  t  Xa  -  1)^u^  +  4v^' I  . 

In  addition,  if  u  £  tf(u  ),  then  the  Rankine-Hugonlot  condition  gives 
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NOTE  1 !  Arrows  on  rarefaction  curves  indicate  the  direction  of  increasing  eigenvalue 

Arrows  on  shook  curves  indicate  the  direction  of  decreasing  shock  speed. 

NOTE  2:  1-  and  2-shocks  in  the  Hugoniot  locus  of  u.  are  indicated  by  dashed  lines 

'  ■ 

supported  by  the  solid  line  for  the  Hugoniot  locus. 


Flerure  Captions 

Figure  1i  Integral  curvea  for  the  eymmetrlc  aysteaa  In  Region  IV>  l.e.,  ayetem  (3) 
with  b  -  0,  a  >  2>  (Family  t  is  denoted  by  the  thicker  linea.) 

Figure  2:  The  Sectors  A2  determined  by  the  ray  6  •  d*  (X^  *  0). 

Figures  3A-F:  Hugonlot  loci  and  shock  types  for  u,  in  representative  positions. 
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Figure  4A:  IRS-composite  wave  in  x,t-8pace 
Figure  48:  2SR-composite  wave  in  x.t-space 
Figure  5:  Typical  curves: 

Rarefaction  (R^(u)  “  luFEj), 

1- Composite  (C,(F)  "■  (EDl ) , 

2- Boundary  ( [“CC*] ) , 

Triple  Shock  ((ABC*)). 

Figure  6A:  The  wave  curves  for  u^^  e  A-j. 
Figure  6B:  The  wave  curves  for  u^^  t  A2. 
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Figure  6b ;  The  Wave 
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We  are  interested  in  classifying  the  solutions  of  Riemann  problems  for  the 
2^2  conservation  laws  which  have  homogeneous  quadratic  flux  functions.  Such 
flux  functions  approximate  an  arbitrary  2^2  system  in  a  neighborhood  of  an 
isolated  point  where  strict  hyperbolicity  fails.  This  problem  was  motivated  by 
Marchesin  and  Paes-Leme  who  discovered  such  a  singularity  in  a  system  of  equa¬ 
tions  arising  in  oil  reservoir  simulation  [15,17].  in  [18],  Schaeffer,  Shearer, 
Marchesin  and  Paes-Leme  solved  the  Riemann  problem  for  this  system  in  a 
neighborhood  of  the  singular  point.  In  [8] ,  Isaacson  and  Temple  outlined  a 
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program  for  classifying  such  singularities  by  means  of  locating  normal  forms 
for  the  equivalence  classes  of  equations  generated  by  linear  changes  in 
dependent  variables.  A  2-parameter  family  of  such  normal  forms  was  found  by 
Plohr.  In  the  important  work  of  Schaeffer  and  Shearer  [17] ,  a  new  normal  form 
was  found  which  reduced  the  classification  of  integral  curves  to  a  theorem 
of  Darboux  on  the  classification  of  umbilic  points  for  homogeneous  cubic 
equations  [3] .  The  integral  curves  fall  into  four  isomorphism  classes,  called 
Regions  I-IV.  In  this  paper  we  give  the  solution  of  the  Riemann  problem  for 
the  systems  in  Region  IV  which  exhibit  up-down  symmetry.  A  presentation  of 
the  solutions  of  the  corresponding  systems  in  Regions  II  and  III  will  follow. 

Our  analysis  uses  a  numerical  determination  of  the  Hugoniot  loci. 
Moreover,  a  new  phenomenon  occurs  in  Regions  Il-IV  that  does  not  occur  in 
Region  I:  lines  exist  on  which  one  of  the  eigenvaJ.ues  has  the  value  of  the 
eigenvalue  at  the  singulcir  point.  Bifurcations  occur  at  these  lines,  and 
the  dynamics  of  the  lines  as  the  parameters  in  the  normal  form  are  veuried 
gives  a  geometric  interpretation  of  the  boundciries  between  Regions  I-III. 

The  program  for  classifying  hyperbolic  singular  points  in  2^2  systems  of 
conservation  laws  is  being  carried  out  jointly  by  the  authors  named  above. 


